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Abstract 

The radiation from the mixed layer into the interior of the ocean of near-inertial 
oscillations excited by a passing storm in the presence of the beta effect is reconsidered 
as an initial- value problem. Making use of the fact that the mixed layer depth is much 
smaller than the total depth of the ocean, the solution is obtained in the limit of an 
ocean that is effectively infinitely deep. For a uniform initial condition, analytical 
results for the velocity, horizontal kinetic energy density and fluxes are obtained. The 
resulting decay of near-inertial mixed layer energy in the presence of the beta effect 
occurs on a timescale similar to that observed. 



1 Introduction 

There is much observational evidence, starting with Webster (1968) and Pollard and Millard 
(1970), that storms can excite near-inertial currents in the mixed layer of the ocean. This 
phenomenon is evident in observations from the Ocean Storms Experiment (D'Asaro et al. 
1995, Levine and Zervakis 1995, Qi et al. 1995). Simple models which treat the mixed layer as 
a solid slab have been quite successful at explaining the process by which wind generates such 
currents (see, e.g., Pollard and Millard (1970), D'Asaro (1985)). A weakness of the model of 
Pollard and Millard (1970) is that it explains the decay of these currents with an arbitrary 
decay constant. Much subsequent work has attempted to determine the detailed character- 
istics of this decay, with possible mechanisms including nonlinear interactions which transfer 
energy to other frequencies (Henyey et al. 1986), turbulent dissipation (Hebert and Mourn 
1993), and the radiation of downward propagating near-inertial oscillations (NIOs) excited 
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by inertial pumping into the interior of the ocean (Gill 1984). The downward radiation of 
NIOs will be the focus of this paper. 

Observations give a timescale for the decay of the energy deposited by the passing storm 
on the order of ten to twenty days (D'Asaro et al. 1995, Levine and Zervakis 1995, Qi et al. 
1995). This timescale stands in contrast with estimates such as that by Gill (1984) that near- 
inertial currents decaying through the downward propagation of NIOs and with a horizontal 
length scale typical of the atmospheric forcing mechanism can remain in the mixed layer for 
longer than a year. To account for this difference, several mechanisms for the enhancement 
of vertical propagation of NIOs have been suggested. D'Asaro (1989) demonstrated that the 
/3-effect causes a reduction of horizontal scales because the meridional wavenumber evolves 
according to I = l — j3t, where l is the initial wavenumber, and / < corresponds to 
southward propagation; this accelerates the rate of inertial pumping of energy out of the 
mixed layer, thereby enhancing the decay. The decay is also enhanced through interaction 
with background geostrophic or quasigeostrophic flow (e.g. Balmforth et al. 1998, Balmforth 
and Young 1999, and van Meurs 1998). 

This paper reconsiders the vertical propagation of near-inertial energy deposited into the 
mixed layer by a storm, in the presence of the /3-effect, using a different approach from that 
of D'Asaro (1989). The analysis uses the formalism of Young and Ben Jelloul (1997) which 
is outlined in Section 2. In Section 3, a simplified model with three main assumptions is 
presented. First, the background flow is assumed to be constant in the zonal direction (i.e. 
independent of longitude with zero vorticity). Second, the buoyancy frequency is taken to be 
small in the mixed layer, and constant in the ocean interior (i.e. beneath the mixed layer). 
Third, it is assumed that the storm has moved very rapidly across the ocean and has created 
a horizontally uniform near-inertial current to the east concentrated within the mixed layer: 
it is the subsequent evolution of this motion that is examined. Section 4 uses the fact that 
the depth of the ocean is very much larger than the mixed layer depth to formulate and solve 
the model for an ocean which is effectively infinitely deep. Section 5 discusses the results 
and suggests directions for further investigation. 

2 The NIO equation 

We consider an ocean of infinite horizontal extent and depth D, with the mixed layer compris- 
ing the region — i7 mix < z < 0, and the rest of the water column occupying — D < z < —H mix . 
The x and y axes are taken to point to the east and north, respectively. The buoyancy fre- 
quency N = N(z) is an arbitrary piecewise continuous function of depth z. 

Young and Ben Jelloul (1997) derive an evolution equation for a complex field A(x, y, z, t) 
which governs leading-order NIO motion in the presence of a steady barotropic background 
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flow and the /3-effect: 

LA t 

where 



^) + i/ V-A + i(^ + lc)LA = 0, 



LA= o_ ffldj; 

dz \N 2 dz , 



(2) 



■0 is the streamfunction for the background flow, ( = V 2 ip is the associated vorticity, and 
the Coriolis parameter is / = / + f3y. Here V is the horizontal gradient, and V 2 = d 2 + d 2 . 
Subscripts denote partial differentiation. The NIO velocity field (w, v,w), buoyancy b, and 
pressure p are given by 

u + iv = e~ ifot LA, 
1 



it; 



V 



--tiN-\A xz -iA yz )e 



+ C.C., 



-fo(A xz - iA yz )e- tfi>t + c.c, 
UA X - iA y )e- ifot + c.c. 



The buoyancy b is related to the density p by 

1 r z 



P = Po 



1 - 
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where p is the reference density at the top of the ocean. The pressure p has been normalized 
by po- 

The boundary conditions are that A z = at z = and z = —D. This ensures that w 
vanishes at the top and bottom of the ocean. Using these boundary conditions, 



/°< 

J-D 



[u + iv) 



0. 



(3) 



Thus barotropic motion is not included in the analysis. However Gill (1984) has shown that 
the barotropic response to a storm is instantaneous and the associated currents are weak. 



3 A Simplified Model 

To simplify the analysis, we assume that A and ip do not vary in the x-direction, and that 
C = 0. The analysis thus neglects the effect of background barotropic vorticity but crucially 
keeps the (3-eSect. The buoyancy frequency profile is taken to be 

N 2 = e 2 N 2 , -H mix <z<0, 
N 2 = N 2 , —D < z < — i/ mix , 
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where e < 1. Finally, the storm is assumed to have produced an initial condition of a 
horizontally uniform near-inertial current to the east concentrated within the mixed layer. 
Instead of approaching this problem by use of an integral operator as in D'Asaro (1989) or 
by projecting onto normal modes (e.g., Gill 1984, Balmforth et al. 1998), the problem will be 
formulated as an initial value problem on a semi-infinite domain corresponding to an ocean 
that is effectively infinitely deep. In order to formulate the problem properly for this limit, 
this section considers an ocean of finite depth. In Section 4 the solution in the limit that the 
depth of the interior is much greater than the mixed layer depth will be found. 

This formulation as a radiation problem which ignores the presence of the ocean bottom 
requires the projection of the initial condition to be spread across all the normal modes. 
This is certainly true for small mixed layer depths in the model of Gill (1984), as shown 
in Table 1 of that paper; also see Table 1 of Zervakis and Levin (1995). For deeper mixed 
layers, this is no longer true since half the initial energy becomes concentrated in the first 
two or three modes. However, as pointed in Section 7 of Gill (1984), the depth of the ocean 
"influences the rate of loss of energy by imposing modulations on the rate, but the average 
rate of loss is not affected very much by depth changes". Hence the results presented here 
should be qualitatively relevant even when the continuum assumption is not valid. 



3.1 Nondimensionalization 

Quantities are nondimensionalized according to 

y = y/Y, z=l + z/H mix , i=ttt, N = N/N Q , 

where 

Typical values (3 = KT 11 nrV 1 , # mix = 100 m, f = 10~ 4 s" 1 , iVo = 10~ 2 s" 1 give Y = 10 5 
m and Q = 10~ 6 s -1 . The relevant timescale is thus Q^ 1 = 11.5 days. Also, the velocity and 
the field A are nondimensionalized by 



(u,v) 



U ' UNgH- 



2 TJ2 

mix 



where U is a characteristic value of the initial velocity. 

The hats are now dropped for ease of notation. With this nondimensionalization, the 
buoyancy frequency profile is 

N 2 = e 2 , < z < 1, 

N 2 = 1, -H = 1 - D/H mbc <z<0, 
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and the NIO equation (|T|), the boundary conditions, and initial condition become 

A zzt + l -N 2 A yy + iyA zz = 0, (4) 

A z = 0, z = -H, z=l, (5) 
A zz = N 2 (u + iv), t = 0. (6) 

The requirement that u and v remain finite imply the jump conditions 

A Z | 2 = 0+ — 'A Z \ Z = Q- , Ayy\ Z = Q+ = Ayy\ Z= Q~, (7) 

where z = + and z = 0~ are the limits as z — > from positive and negative z values, 
respectively. 

This nondimensionalization allows some immediate conclusions to be drawn about the 
propagation of NIO energy downwards. Most importantly, if H mix increases, then the 
timescale decreases. Thus, assuming that the storm causes a uniform near-inertial 
current throughout the whole mixed layer, energy transfer will be faster for a deeper mixed 
layer. This confirms the results of Gill (1984), which associated the more efficient transfer 
with a larger projection of the initial velocity profile on the first vertical mode. 

3.2 Boundary Condition at the Base of the Mixed Layer 

Expanding A(y, z, t) = A Q (y, z, i) + e 2 A 2 (y, z, t) + C(e 4 ) for < z < 1, (g) becomes at O(e ) 

A 0zzt + iyA 0zz = 0. 

Integrating this subject to the boundary condition that A z and thus A 0z vanishes at z — 1 
implies that A Q is independent of z. At 0(e 2 ), 

A 2zz t + iyA 2zz + -A 0yy = 0, (8) 

which may be integrated subject to the boundary condition that A 2z vanishes at z = 1 to 
give 

A 2zt + iyA 2z + -A Qyy (z - 1) = 0. 
Evaluating at z = + and using A yy = A 0yy + 0(e 2 ) and A z = e 2 A 2z + C(e 4 ), 

ie 2 

A zt + iyA z - —Ayy = 0(e 4 ), z = + . 

Finally, applying (|7|) gives the upper boundary condition for the NIO field in the ocean 
interior to leading order in e: 

A zt + %yA z - l -A yy = z = IT. (9) 
5 



Results obtained in the ocean interior using (^) are in fact leading-order solutions. We 
shall continue to use the notation A, even though it is really the leading-order term in the 
expansion. 

3.3 Initial Condition 

Suppose that in a short time compared with the NIO wave propagation time, the passing 
storm induces near-inertial currents in the mixed layer with a horizontal scale that is much 
larger than the one under consideration, and which can hence be taken to be uniform. For 
simplicity, the initial velocity (consistent with equation (|3])) is assumed to be piecewise 
constant with depth: 

(u,v) = (1,0) < z < 1, 

= (-tf-\0), -H<z<0. 

The weak flow in the ocean interior is necessary to ensure that the flow has no barotropic 
component. Integrating equation with respect to z and using the boundary conditions 
(|) gives at t = 

A z = e 2 (z-l), 0<z<l, (10) 
A z = -(z + H)/H, -H<z<0. (11) 



4 Solution for an Infinitely Deep Ocean 

The total depth of the ocean is typically on the order of a hundred times the depth of the 
mixed layer. Thus, the limit of infinite depth is considered. The initial condition is taken 
to be equation (|ll|) with H — > oo. The boundary condition for z — > — oo is taken to be 
A zz — > 0, corresponding to the near-inertial velocities vanishing at infinite depth. Of course, 
this limit excludes the possibility of reflections off the bottom of the ocean which may be 
important. Finally, the boundary condition for z = 0~ given by equation @ is used. Hence 
the problem to be solved for the semi-infinite domain z < becomes 

Kzt + -A yy + iyA zz = 0, z < 0, 

A zt + iyA z - -A yy = 0, z = 0~ , 

A zz ^0, z -> -oo, 

A z = -1, t = 0. 
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4.1 NIO velocity field 



These equations may be solved using Laplace transforms. Here we present only the major 
results; further details are given in Moehlis (1999). We make the transformations A(y, z, t) = 
e~ iyt B(z,T), T = t 3 /3, and a = (1 + i)/2 and define the Laplace transform of B by 



Then 



b(z,p) = C[B] = / B(z,T)e- pT dT. 
Jo 



K z ,p) = — —^- A — ex P 



(12) 



(13) 



This Laplace transform and its derivatives with respect to z must be inverted numerically 
for the ocean interior (z < 0). For the top of the ocean interior (z = 0~) however, they may 
be obtained in closed form. For example, 



A ez (y,0 ,t) = e 



-iyt 



(14) 



We now consider the back-rotated velocity A zz = e 4 -' {u + iv), which filters out purely 
inertial motion at frequency fo. Back-rotated velocities may be represented by hodographs 
which show the vector (Re(A zz ), Im(A^)) as curves parametrized by time. For f > 0, 
if these curves are traced out in a clockwise (counterclockwise) fashion, the corresponding 
motion has frequency larger (smaller) than f . Figure [I] shows the back- rotated velocity 
at different locations. A common characteristic is that the magnitude of the back-rotated 
velocity starts at zero, reaches a peak value shortly after the storm, then decays away. The 
depth dependence of the back-rotated velocity is seen by comparing Figure [3] (a) and (b), 
where both have y = and thus the same value of the Coriolis parameter /. Qualitatively the 
results are the same, but closer to the mixed layer the direction change of the back-rotated 
velocity becomes slower, meaning that the frequency is closer to fo. An idea of the latitudinal 
dependence is seen by comparing Figure [1] (a,c,d): at y = 1 the hodograph is traced out in a 
clockwise fashion as for y = 0, but at y = —2 it is traced out in a counterclockwise fashion. 



4.2 Kinetic energy density and fluxes 

The horizontal kinetic energy (HKE) per unit area contained within the mixed layer is 



dz 



A ? 



N 2 



= dz 
Jo 



A, 



I dz \A 2z 
Jo 



Expanding B(z,T) = B (z,T) + e 2 B 2 (z,T) + C(e 4 ) in the mixed layer, (||) may be used to 
show that 



pb 2 zz - B 2zz (z, 0) - -b = 0, 



(15) 
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Figure 1: Back-rotated velocity for (a) z = — 1, y — 0, (b) z — —0.5, y — 0, (c) z — — 1, 
y = 1, and (d) z — — 1, y — —2. The diamonds are drawn at t — 0, 5, 10, 15, 20. 
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Figure 2: Horizontal kinetic energy per unit volume (HKE) in the mixed layer, cml, for 
linear and logarithmitc axes. The solid line shows the exact result and the dashed line the 
asymptotic result. 

where b 2 = C\B^ and b = C[B ]. The initial condition within the mixed layer is B 2zz (z, 0) = 
1. Now A is continuous across z = 0, and B is independent of z (see Section 3.2). Hence 



1 



>1zz 



p 2ap ^fp + a ' 



which may be inverted to give 



A 2zz {y,t) = e-^e^erfc . 



(16) 



Therefore the HKE within the mixed layer is 



erfc 



f l + i 



f 3/2 



The time dependence of cml is shown in Figure [| Asymptotic results from Abramowitz 
and Stegun (1972) for the complementary error function imply that 



1 - 
6 



t — > OO. 



t < 1, 
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Since the energy which leaves the mixed layer enters the interior of the ocean, this implies 
that for short times the energy in the interior increases like t 3//2 . This does not contradict the 
result from D'Asaro (1989) that for short times the thermocline energy grows like t 6 . That 
result assumes that the wind persists to generate a constant inertially oscillating velocity, 
and that there is no propagating inertial motion. Here, the wind has an instantaneous effect, 
causing an initial horizontally uniform inertial current, and propagating inertial motion is 
included fully 

Another quantity of interest is the flux of HKE. Using (Q) and its complex conjugate 
gives 



d d 

s HKE =s 



.4 



zz 



i d , . ... _ , , i d 



2PP ~dy^ AzzA * y ~ A * zzA ^ + 2iV2 'dz <(A * yzAy ~ A y zA v">- ( 17 ) 



N 2 

Assuming A ZZ A* — A* zz A y vanishes for \y\ — > oo and using equation 



d 



/—d poo poo poo poo 

dz dx dy \A ZZ \ 2 = / / d) dx dy, 

-H J —oo J —oo J —oo J —oo 



dt J-H 

where 

i 

2' 

gives the flux of HKE from the region z > —d to the region z < —d. For this model, we 



F E (y,t;d) = -(A* yz A y - A yz A* y )\ z= _ d (19) 



consider the flux per unit area. Integrating fll8D with respect to time shows that the quantity 



E(t;d) = f F E dt 
Jo 



gives the total amount of HKE which has penetrated into the region z < —d. Note that 
E(t; d) — > 1 corresponds to all the energy originally in the mixed layer having reached 
depths below z = —d. Results for F^it] d) and E(t; d) obtained by numerically inverting the 
appropriate Laplace transforms are shown in Figure |3|. Fe peaks at the nondimensionalized 
time t w 0.62; for the typical values quoted in Section 3.1, this corresponds to about a week 
after the storm. From Figure |3](b) and using the fact that whatever energy flows through 
z = 0~ must have initally been in the mixed layer, we see that by t — 1 (about 11.5 days 
after the storm) nearly half of the energy associated with horizontal NIO currents caused 
by the storm has left the mixed layer; however, only about 38% of the total energy has 
penetrated below z = — 1. Byt = 2 (about 23 days after the storm), 82% of the total energy 
has left the mixed layer, but only 58% has penetrated below z = — 1. Thus, at t = 2 nearly a 
quarter of the total energy is contained in the distance H m i x immediately beneath the mixed 
layer. This is reminiscent of the accumulation of NIO energy below the mixed layer seen in 
Balmforth, Llewellyn Smith and Young (1998). This model thus gives reasonable estimates 
for the timescale for which the decay of NIO energy occurs: for example, D'Asaro et al. 
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Figure 3: (a) Fg(t;c?) and (b) E(t;d) for different depths d below the base of the mixed 
layer. These show instantaneous and time-integrated fluxes of HKE. 



(1995) found that the mixed layer inertial energy was reduced to background levels by 21 
days after the storm. 

Figure [| shows the vertical dependence of the HKE and Fe at different times. As time 
increases the instantaneous distribution of HKE becomes more sharply peaked near the base 
of the mixed layer, but remains bounded (asymptotically approaching unity) because of 
energy conservation. 



4.3 Large- time behavior 

The asymptotic behavior of near-inertial properties may be derived using the method of 
steepest descents (see Moehlis 1999 for details). This shows that in the limit of large £ = 
z 2 / 3 t, and along the "rays" z = — ?7o^ 3 /3, 

2,2 2 e, 2 Vo 

(1+^)^' 7r(l+IJ§)f 

A useful way to represent the asymptotic results is to write r/o in terms of z and t and then 
draw contour plots of quantities of physical interest in the (z, t) plane: this is shown in 
Figure |5|. In the asymptotic limit for large £, with z constant, u 2 + v 2 and Fe decrease as 
time increases. Note that £ is large for sufficiently large z and/or t. 

Finally, Moehlis (1999) also obtained results for the vertical shear u 2 + v 2 . To leading 
order in e, the vertical shear within the mixed layer is zero. The results for vertical shear for 
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Figure 4: Vertical profiles of (a) u 2 + v 2 and (b) -Fe(£, \z\) at y = for different times showing 
the decay of energy from the mixed layer (0 < z < 1) and resultant behavior in the interior 
[z < 0). Note the different vertical scales. 

the interior of the ocean lack physical realism because the model allows the shear to grow 
forever as a consequence of the initial infinite shear due to the discontinuity in the initial 
velocity profile. 



5 Conclusion 

A simplified model has been developed to examine the decay due to the /3-effect of near- 
inertial currents excited in the mixed layer by a passing storm. This decay occurs due 
to the radiation of downward propagating NIOs into the interior of the ocean. The main 
assumptions of the model are that the background flow does not vary in the longitudinal 
direction and has no associated vorticity, that the ocean has a simple (piecewise constant) 
buoyancy frequency profile, and that the storm has moved very quickly over the ocean causing 
a horizontally uniform near-inertial current concentrated in the mixed layer. The /9-effect is 
included in the analysis and is responsible for the radiation of NIOs. Because the depth of 
the mixed layer is much smaller than the total depth of the ocean, the problem is formulated 
in the limit of an effectively infinitely deep ocean; the resultant initial value problem is solved 
by Laplace transforms. Analytical results are given for the horizontal kinetic energy density 
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2 2.5 3 3.5 4 4.5 5 

Figure 5: Contour plots of the asymptotic results for (a) u 2 + v 2 and (b) Fe- Darker shading 
corresponds to smaller values. 

in the mixed layer, and results from the numerical inversion of the appropriate Laplace 
transforms are given for horizontal kinetic energy, energy flux, and back-rotated velocity. 
The asymptotic behavior is also investigated. 

Although this simplified model cannot be expected to capture the full complexity of the 
aftermath of a storm passing the ocean, it does capture much of the observed behavior. 
Most importantly, in the presence of the /^-effect the decay of near-inertial mixed layer 
energy is found to occur on the appropriate timescale (approximately twenty days), which 
confirms the analysis of D'Asaro (1989) and observations by D'Asaro et al. (1995), Levine 
and Zervakis (1995), and Qi et al. (1995). The main advantage of the approach described in 
this paper is that many aspects of the decay in the mixed layer are analytically obtained for 
all times, unlike D'Asaro (1989) which predicts the timescale for the decay in a short time 
limit or estimates it in terms of the time it takes normal modes to become out of phase (cf. 
Gill 1984). Extensions to a more realistic ocean and storm would involve including a more 
realistic buoyancy frequency profile (for example, the profile used by Gill 1984), considering 
the effect of different initial velocities (including both horizontal and vertical structure), and 
considering the effect of background flow. The study of all of these could use the same 
formalism of Young and Ben Jelloul (1997) and an approach similar to that presented here. 
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